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Abstract 

The conical self-similar vortex solution of Long (1961) is reconsidered, with a view toward understand- 
ing what, if any, relationship exists between Long's solution and the more-recent similarity solutions of 
Mayer and Powell (1992), which are a rotational- flow analogue of the Falkner-Skan boundary-layer flows, 
describing a self-similar axisymmetric vortex embedded in an external stream whose axial velocity varies 
as a power law in the axial (z) coordinate, with <j> = r/z n being the radial similarity coordinate and n 
the core growth rate parameter. We show that, when certain ostensible differences in the formulations 
and radial scalings are properly accounted for, the Long and Mayer-Powell flows in fact satisfy the same 
system of coupled ordinary differential equations, subject to different kinds of outer-boundary conditions, 
and with Long's equations a special case corresponding to conical vortex core growth, n — 1 with outer 
axial velocity field decelerating in a z~ x fashion, which implies a severe adverse pressure gradient. For 
pressure gradients this adverse Mayer and Powell were unable to find any leading-edge-type vortex flow 
solutions which satisfy a basic physicality criterion based on monotonicity of the total-pressure profile of 
the flow, and it is shown that Long's solutions also violate this criterion, in an extreme fashion. Despite 
their apparent nonphysicality, the fact that Long's solutions fit into a more general similarity framework 
means that nonconical analogues of these flows should exist. The far-field asymptotics of these general- 
ized solutions are derived and used as the basis for a hybrid spectral-numerical solution of the generalized 
similarity equations, which reveal the existence of solutions for more modestly adverse pressure gradients 
than those in Long's case, and which do satisfy the above physicality criterion. 

1 Outline 

A brief outline of the paper is as follows: We first review the Hall and Mayer-Powell leading-edge-type and 
Long "tornado-like" similarity solutions. We then show how both of the different similarity procedures used 
for these disparate flow types can be reduced to a single unified framework accommodating a continuum of 
growth-rate-parameterized solutions, in which what we dub "Long- type" solutions satisfy the same similarity 
equations derived by Mayer and Powell but with a far-radial-field at infinity. The particular flows described 
by Long are a special conical- vortex-growth case of the general solutions. We next derive asymptotic series 
expansions describing the behavior of the generalized Long-type flows in the far radial field, and describe 
a novel hybrid spectral-numerical procedure for solving the corresponding full similarity equations in the 
inner viscous core region which automatically includes the matching to the appropriate far-field asymptotic 
solution, by way of treating the unknown far-field series coefficient as an eigenvalue of the problem. The 
numerical results show that for less-adverse pressure gradients than those in Long's case (and thus having 
viscous core growth rates less than conical), these full- field solutions do satisfy the above physicality criterion 
based on total-pressure. The asymptotic analysis of the far-field behavior of the generalized solutions also 
shows that there is a close correlation between the onset of nonphysicality according to the total-pressure 
criterion and singular behavior in the coefficients of the asymptotic series. The singularity is only apparent 
when one considers a continuum of solutions with different values of the growth-rate parameter, and thus 
would not be apparent if a single similarity law (e.g. conical growth, as in Long's flows) is considered. 



2 Introduction 

In this study, we consider an idealized model of axisymmetric vortices with slender cores in an incompressible 
fluid having infinite extent. Under these assumptions, two kinds of vortex flows are frequently studied. The 
first kind is a leading-edge-type vortex with nonzero axial velocity excess over that in a "far field" taken to 
be sufficiently outside the immediate vicinity of the viscous core but not extending to infinity; herein we will 
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refer to such flows as "Hall-type vortices" for simplicity. The second kind consists of a viscous swirling core 
embedded in a potential vortex, i.e. having circulation r = vr tending to a constant in the far field, and all 
velocity components decaying to zero far from the axis; in the context of the present study we shall refer to 
these type of solution of the governing equations as "Long-type vortices." 



[Hall (1961) proposed a similarity-solution model for vortices formed by the rolling up of the shear layer 
produced at the leading edge of a lifting delta wing, which divided the vortex core into two regions: an 
inviscid rotational conically self-similar outer core, with viscous effects appreciable only in a nonconical 
inner subcore. This approach yields an exact solution for the conical outer region with velocity components 
(it = radial, v = azimuthal, w = axial) of the form 
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where 4> = r / z i a = \/l + 2V^/W^ — 1 and 4> e is an outer-boundary-condition location whose value can 
be chosen - the only obvious constraint is that it should be sufficiently far from the rotational axis that 
viscous effects are negligible there, - a requirement which can be checked a posteriori - and V e and W e 
are the azimuthal and axial velocity values (also user-specified) at this outer edge location. For the inner 
viscous core, Hall was able to find a solution by making scaling assumptions of the usual boundary-layer 
type, subject to two key additional assumptions: (i) the axial velocity within the viscous subcore is nearly 
constant (which is not justified physically, based on the large axial-velocity variations seen in such flows in 
experimental studies such as the contemporaneous one of |Earnshaw( 196 1 )] , and (ii) the radial velocity has 
a far-field near-linear velocity profile which matches the inviscid outer core solution. With these additional 
assumptions, Hall was able to reduce the governing equations to an analytically tractable form. The inner- 
core solutions obtained this way are in fairly good agreement with Earnshaw's experimental data in the outer 
part of the viscous subcore region, but give predictions for the inner core (the region of greatest interest) which 
only agree with the experimental data in a rough qualitative sense. Hall's approach additionally made it 
difficult to examine the solution behavior in the limit of vanishing viscosity, and utilized an unsatisfactory ad 
hoc matching procedure between the inner and outer solutions. Soon thereafter, Stewartson fc Hall(1963)| 



improved the inner subcore solution via a proper asymptotic matching approach; the actual results, however, 
were little different from Hall's original work. Fundamentally, due to the nature of the above outer core flow, 
the solutions can only be extended a finite distance from the vortex axis irrespective of which matching 
technique one uses. This is perhaps not terribly troubling from the standpoint of understanding leading- 
edge- type vortex cores since such flows naturally have a limited radial extent (since the vorticity is introduced 
at a finite radial distance by the shear layer rolling up from the leading edge), but it does beg the question of 
whether there is any kind of flow model of this type which leads to solutions which can be extended radially 
to infinity, as those for wake-like and potential vortices can. 

Mayer & Powell(1992)j generalized Hall's approach by allowing for variation of axial velocity in the 



far field of power-law form W(z) ~ z m . They made the added simplification of using a single similarity 
variable for the entire core region in lieu of the separate "inner" and "outer" core approach used by Hall. A 
straightforward similarity analysis showed that the viscous core region has an axial growth rate described 
by a parameter n which appears by way of the radial similarity variable <fi = r/(ez n ) and is coupled to the 
power-law parameter that describes the axial variation of the far field via the simple relation n — 1 ~ 2 m . (The 
power-law scalings for these flows are in fact identical to those of the well-known Falkner-Skan boundary- 
layer flows, which analogously generalize the famous solution of Blasius (1908) for the viscous boundary layer 
formed by flow over a flat-plate with constant far-field conditions.) They thus obtained a set of similarity 
equations describing a generalized family of Hall-type flows. 

The equations for m = 0, i.e. with no axial velocity or pressure gradient and with nonconical self- 
similarity, admit an analytical solution in the outer core region which is identical to Hall's outer solution ([1}, 
but with the radial velocity u multiplied by the growth rate parameter n. Owing to the single similarity 
variable, the full viscous solution can be obtained numerically to any desired accuracy without requiring a 
tedious matching process between the inner and outer core regions, and the m = solutions agree quite well 



with the high-quality experimental data of Verhaagen & Ransbeeck(1990) for incompressible flow over a 



delta wing. For m decreasing from 0, which corresponds to an increasingly adverse axial pressure gradient 
of the outer flow, the axial velocity profile of the solutions shows an increasingly pronounced retardation in 



2 



the viscous region near the vortex axis, flow behavior which is reminiscent (within the obvious limitations 
of the similarity approach) of the large-scale flow disruptions seen in the phenomenon of vortex breakdown 
in real- world flows (cf. Sarpkaya(1971)| ). 



Note that in both the original Hall-type formulation and its generalizations, since one is specifying 
boundary conditions at a finite distance from the vortex axis, unless one manages to eliminate the pressure 
from the formulation, this includes specifying a pressure value at the outer boundary. For m — the pressure 
appears in the similarity equations only via its derivative, hence this outer-boundary value is irrelevant in the 
sense that neither the velocity profiles nor the degree of pressure change between the outer flow and the vortex 
core depend on the chosen boundary value. However, for nonzero values of m an undifferentiated pressure 
term appears in the similarity equations, as well, which means that the chosen value of edge pressure has a 
deterministic effect on the solution properties. Mayer and Powell dealt with this issue by showing that there 
is a particular choice of edge pressure such that the resulting distributions of total pressure p + \(v 2 + w 2 ) 
(in terms of the similarity scalings, the radial velocity u is negligible and drops out of this expression) 
are asymptotically constant and decrease to a minimum on the axis, which is physically realistic based on 
behavior of real flows. These finite-boundary-distance- related issues are absent in flow models in radially 
unbounded domains such as the Long- type flows we shall describe next, even though (as we will also show) 
the governing similarity equations are the same. 

Dating back to virtually the same time period as the leading-edge-vortex work of Hall, the family 
of conically self-similar rotational solutions of the steady incompressible Navier-Stokes equations discov- 
ered by Long(1961)| and collectively referred to as "Long's vortex" is commonly used as an idealized 



model of tornado-like flows ( [Burggraf fe Foster(1977)| Shtern fc Hussain(1993)| ) and as a reference flow 



for studies of swirling-flow stability (poster fc Duck(1982)| |Foster fc Smith(f 989 j[ [Foster fc Jacqmin(1992)| 



Khorrami & Trivedi(f 994) Ardalan, Draper fc Foster(1995)| ). Because of this wide interest, one would as- 



sume that the basic flows themselves have been subjected to close scrutiny, found to be physically reasonable 
according to some rational set of criteria, and the details of their scalings understood and placed in some 
clear relation to other boundary-layer-like flows, of both the swirling and non-swirling variety. We find that 
this is not the case. Some important points left unaddressed in this regard in the extant literature are: 

1. In his derivation, Long used a "boundary- layer-like" (his words) scaling for the radial similarity variable 
of y ~ r/(ez,z), where €l = Re -1 and the Reynolds number is based on outer-flow circulation. Based 
on this definition, the radius of the core grows linearly with axial distance downstream of a singular 
point at the apex of the vortex, and is inversely proportional to Reynolds number. Compared to the 
0(Re~ 5 ) thickness growth rate of other well-known boundary-layer flows, the scaling used by Long is 
unusual. (In fact it is precisely this odd scaling that first led us to reinvestigate Long's formulation.) 

2. Perhaps the most prominent feature of Long's solutions is their nonuniqueness with respect to a charac- 
teristic parameter </, the so-called "flow force," equivalent to the total flux of axial momentum through 
any plane perpendicular to the vortex axis. For any given value of J in a particular range, there exist 
two sets of solutions of the similarity equations, which Foster and Duck (and subsequent authors) la- 
beled type I and type II flows, and which are mainly differentiated by their axial velocity components: 
the type I flows exhibit regions of reversed axial flow, whereas the type II flows do not. We shall show 
that the solutions can in fact be alternately and uniquely parameterized by their axis values of axial 
velocity, and thus that the apparent nonuniqueness is strictly a result of Long's particular choice of 
parameterization. However, from a physical standpoint, no plausible physical mechanism has been 
advanced for the reversed flows. 



At large values of the flow force, Foster fc Smith(I989)] found flow variations on the reversed-flow 



(type II) solution branch to become intensely concentrated in an asymptotically narrow annulus about 
r = An 2 J / (pT 2 ) where axial and azimuthal velocities are large, a flow pattern they dub a "ring-jet." 
Peculiar as such solutions are, no questions were raised about their physical plausibility. 

In Long's solutions, the static pressure drops monotonically from a constant value in the far field to 
a minimum on the axis, as expected on physical grounds for a swirling flow. However, to the best of 
our knowledge, no total-pressure plot for any of these flows has previously appeared in the published 
literature. If one neglects the temperature change caused by viscous dissipation (as is justified in the 
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incompressible- flow limit), the total pressure represents the variation of total flow enthalpy. As there 
is no energy source embedded in the flow, one expects the total pressure to reach a minimum on the 
vortex axis (where the cumulative dissipation is greatest) and to increase toward a constant value in 
the far field. This is a minimum requirement for such a flow to be physically plausible. We find that 
none of Long's solutions satisfy this simple total-pressure-based physicality criterion, and that it is 
violated particularly badly by the type I solutions, those with a strictly jet-like axial flow profile 0. 

Due to these unresolved issues with respect to the Long's flows, we decided to undertake a thorough 
review of the basic formulation and solutions for the Long-type vortex flows, not only to to be able to say 
something about their physical plausibility (or lack thereof), but also in an effort to try to place them in 
some kind of context among other well-known vortical-flow similarity models. The gist of our findings, and 
an outline of the paper, is as follows: 

In Sj3] we show that the similarity equations derived by Long are in fact the same as those obtained 
via more-orthodox boundary layer scalings. A corollary of this analysis is that, when properly scaled, 
the similarity laws and equations governing Long's vortex arc precisely the same as those derived by 
Mayer fc Powell(1992)| in their investigation of leading-edge-type swirling flows, but subject to a different 
set of far-field boundary conditions. Placing Long's flows in this context also makes clear why reversed-flow 
solutions might be expected to exist: conical self-similarity, as assumed a priori by Long, is only possible if 
the outer axial flow decelerates in a z~ l fashion, which implies a highly adverse axial pressure gradient, in 
fact one which is much more severe than the most-adverse case for which the aforementioned authors were 
able to find physically reasonable leading-edge-type solutions, and similar nonphysicality is evinced by Long's 
solutions as well, once one examines their total pressure profiles. Fortunately, all is not lost. Since the Long's 
flows are solutions of a special case (having m = —1) of the Mayer-Powell similarity equations, there exists 
the possibility that there are generalized Long's flows of nonconical form, most especially ones at less-adverse 
axial pressure gradients (m > — 1) and hopefully having more-reasonable total pressure distributions. The 
derivation of these generalized (nonconical) analogues of Long's similarity equations concludes the section. 

Prior to attempting a numerical solution of the generalized similarity equations, it is worth examining 
(and in fact necessary in the context of the spectral-numerical approach used here) the asymptotic behavior 
of solutions in the far field, and this is the subject of SJU It i s shown that solutions of the equations for 
which all velocities decay at radial infinity are only possible if the pressure gradient is adverse (m < 0) in 
terms of the Falkner-Skan-type similarity parameter governing the flow. General asymptotic series suitable 
for describing the dependent variables (stream function and circulation) are developed, which reduce to the 
well-known far-field series governing Long's solutions in the special case of conical self-similarity, but which 
also prove suitable for nonconical flows. 

Information about the viscous core region appears in the asymptotic far-field series as an undetermined 
coefficient in terms of which the coefficients of all the higher-order terms can be expressed, and which in 
a classical matched-asymptotic-expansion approach would be successively approximated by matching to a 
functional hierarchy which approximates the viscous inner solution. We adopt a rather different approach to 
the matching by treating the undetermined coefficient as an eigenvalue to be found in a numerical solution 
spanning the entire flow field, which is constructed (using the preceding far-field series expansion) to have 
the proper asymptotic behavior in the far radial field, assuming the "matching constant" yielded by the 
numerical solution procedure converges as the numerical resolution is increased. This procedure is described 

1 A reviewer has objected to this criterion as being unjustified under the assumptions used in the flow model, and indeed, 
based on strictly icnompressible flow, it is, as there is no energy equation, no second law of thermodynamics, and hence no 
dissipative "losses" which should be reflected in added constraints on the solutions. But in fact theoretical fluid mechanics 
abounds with such apparent paradoxes and solution selections "unjustified by the model". A classical example is the idealized 2- 
D model of the circular spinning cylinder used by such luminaries as Rayleigh to illustrate crucial concepts related to circulation 
and lift as embodied by the Magnus effect. If the cylinder is truly of a perfectly circular cross-section with a perfectly smooth 
surface and the fluid is truly inviscid there is of course no way for the fluid to "know" whether the cylinder is spinning or not. 
The "unjustified" thing done is to invoke the viscous vorticity transfer known to occur in real fluids as a result of the no-slip 
condition. Does this mean we should ignore Rayleigh's work and conclusions? 

We feel it important, even in the context of a highly idealized athematical model of the physics of interest, to emphasize 
those solutions of the model equations which appear to be of greatest relevance to phenomena observed in the real world. In 
the present case, we use knowledge of total-pressure profiles seen for similar flows in actual experiments to restrict our focus. 
The generalized similarity formulation certainly admits of a similarly wide variety of "heretofore unseen in the laboratory or 
in nature" solutions as the does the original conical Long's formulation. Readers interested in studying these wilder regions of 
the solution space are welcome to do so to their hearts' content. 
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in fJSl as is the high-accuracy iterative spectral method used to solve the resulting system of coupled nonlinear 
ordinary differential equations containing the necessary eigenvalue term. 

The character of the generalized flows is the subject of SjU We begin by reviewing the properties of the 
classical conical case considered by Long in £16.11 As previously mentioned, many aspects of these solutions 
are well-studied, and we use our numerical results and the accompanying plots mainly to illustrate key 
qualitative aspects of the solutions. The discussion concludes by presenting the total-pressure profiles for a 
flow which typifies the Type I solution branch to demonstrate the likely nonphysicality of the conically-similar 
solutions. 

The characteristics of the generalized Long-type flows are the subject of M6.21 An interesting property 
of the generalized similarity equations is that the form of the higher-order coefficients of the corresponding 
asymptotic series derived in SjJ] permits circumstances under which these coefficients blow up due to a 
vanishing denominator. This can only happen if the flow in question has a growth rate which is less than 
conical, and some examples where this singular behavior actually does occur are detailed. Perhaps tellingly, 
the locations in the parameter space where such behavior occurs correlate closely with the onset of near-axis 
overshoots of the total pressure, so this "coefficients crisis" may serve as a useful proxy for the onset of 
nonphysicality as based on this criterion. In ^6.31 we conclude by explicitly demonstrating how one may 
start with a generalized Long-type solution, cut it off at finite radius, and smoothly transition to a Hall-type 
solution at the same flow parameters by varying the outer-edge pressure value, assuming one is in a region 
of the parameter space (defined by growth-rate parameter n and core axial velocity) for which both types of 
solutions exist. 



3 Governing Equations and Similarity Formulation 



The governing equations are the usual time- invariant incompressible axisymmetric Navier-Stokes equations. 
The only salient difference in this regard between the Hall and Long-type formulations is that in the former 
one typically uses a characteristic axial velocity with which to nondimensionalize the flow variables and 
on which to base the Reynolds number, whereas in the latter one uses the asymptotic value of far-field 
circulation for this purpose. In his formulation Long began with the second (circulation-based) definition of 
the Reynolds number and used that to define a small parameter and radial similarity variable y, with y 
assumed to be of order unity in the viscous vortex core: 



e L := Re" 
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along with a nondimensionalized stream function f(y) and circulation g(y), also taken to be of order unity. 
Long claimed the above definition of small parameter and scaling of the radial similarity variable to be 
"boundary-layer- like," but in a standard boundary layer treatment one instead uses e = Re -1 ^ 2 . However, 
it is easily verified that Long's eventual similarity equations in fact express the same balance of terms 
which results from a more conventional boundary layer approach. In fact, the potential confusion here 
arises because Long's scalings in effect correspond to the use of a characteristic radial (i.e. transverse) 
velocity in the nondimensionalization procedure and definition of the Reynolds number. (Whether the 
resulting flow solutions actual exhibit a radial velocity distribution which allows for a suitable reference 
radial velocity value to be defined is thus rendered moot.) With the further expedient of omission of 
the superfluous v2 factor introduced by Long into the denominator of the radial similarity variable, the 
resulting equations (still in terms of nondimensionalized velocities and pressure) are identical to the conical- 
form special case n — —m = 1 of those obtained by Mayer fc Powell(1992)| based on the radial similarity 
coordinate <p — r/(ez n ): 
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where we leave things in terms of both the velocity-field axial power-law parameter m and the viscous 
boundary-layer growth-rate parameter n for simplicity, even though these parameters are related via m = 
1 — 2n, the same relation as in the classical Falkner-Skan boundary layer flows. 

For a leading-edge-type flow, where the spiral roll-up of the shear layer creates a radially finite vorticity 
field, suitable boundary conditions are 

u = v = w' = at <f> = 0; v = v e , w = l,p = p e = — (1 + v 2 ) at <j> = <j) e . (4) 

The choice of "outer edge" location <fi e and edge swirl v e determine the strength of the resulting vortex. The 
edge pressure boundary condition (which is non-arbitrary for m ^ because in that case the pressure is 
nontrivially coupled to the velocities by way of its appearance in the axial momentum similarity equation) 
is obtained from assuming the flow to be irrotational in the far field {<j) e must be taken sufficiently large to 
ensure that this approximation is reasonable) and applying Bernoulli's law, with the total pressure taken as 
zero in the far field. The reason an analogous issue does not not arise for the Falkner-Skan boundary-layer 
flows is that in the latter, the pressure may be assumed constant in the transverse direction, whereas in 
a swirling flow this is not the case. Of course for m 7^ 0, solutions of the above equations may exist for 
other values of p e than that used above which differ nontrivially (i.e. in more than just an arbitrary additive 
constant) from those using p e = |(1 + vf), but as described by Mayer and Powell, none of them appear to 
satisfy the basic criterion of physicality we use here, namely that the total pressure po — P + \{v 2 + w 2 ) 
tends monotonically toward a constant as cj> — > p e . 

For convenience, we will refer to solutions of the similarity equations subject to boundary conditions 
of the form ((4]) as Hall-type solutions. A second class of solutions, and the main topic of this paper, is 
represented by solutions which extend to infinity (when such exist), and for which all flow variables decay 
to zero in the radial far field, i.e. subject to outer boundary conditions 

it, v, w,p — > as 4> — > oo. (5) 

Solutions satisfying these latter boundary conditions will henceforth be referred to as Long-type solutions. 
Since all velocity components decay to zero for this latter class of flows, one cannot adjust the strength of 
the vortex in the same way on does for Hall-type flows, by fixing nonzero values of v and w at the outer 
boundary; rather, one must do something at the vortex axis, for instance setting the axial velocity there by 
replacing the inner boundary condition w'(0) — with w(0) — w ax is- (Alternatively, one might specify v' on 
the vortex axis, in effect specifying the vortex strength by way of the axial component of vorticity there.) 



3.1 On the parameterization of solutions 

The family of flows first studied by Long have historically been parameterized via a quantity called the "flow 
force," defined as 

J ■= = J ip ~ Poo + -pw 2 )rdr (6) 

and which is essentially the nondimensionalized axial momentum flux of the flow through a plane of constant 
z. The proportional parameter M is the one most often used to parameterize solutions; for instance Foster 
and Duck found no solutions for M < 3.75, and two distinct solutions for each value of M greater than this 
critical value (cf . figure 3 in their paper) . Foster and Duck refer to solutions on the upper and lower part of 
the flow-force curve as Type I and Type II flows, respectively; the Type I flows are distinguished by their 
larger axial velocity maxima, whereas the Type II flows (except for a very small portion of the parameter 
space near the minimum M for which solutions exist) exhibit regions of reversed axial flow in their cores. 
The Type I flows also exhibit axial velocity deficits in their cores for small M, but the deficits become less 
pronounced with increasing flow-force, and for M > 4.71 (as determined by the above authors), the axial 
velocity profiles become strictly monotone decreasing functions of the radial coordinate. The existence of two 
solutions for each value of M greater than the critical one has been accepted as an indication of nonuniqucness 
of solutions, but this interpretation is a parameterization-dependent one: the same aforementioned figure 
used by Foster and Duck - which shows a solution curve in a Cartesian coordinate system having M as 
its abscissa and the core axial velocity w(0) as ordinate - demonstrates that solutions are in fact unique if 
instead parameterized by w(0), and this is the parameterization we shall use here. 
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3.2 Stream-function/circulation formulation 



To show explicitly that equations (I3al -d) contain Long's equations as a special case, and also to simplify the 
far-field asymptotic analysis of the generalized Long-type flows, we now reformulate the similarity equations 
in terms of a stream function /, chosen to automatically satisfy the continuity equation (|3al) . and a circulation 
function g, defined as follows: 

£l f 9 I 

u:=nf-~, W ^lj>' 

For the case of conical flow (n = 1) these definitions are equivalent to Long's, except that Long multiplies 
the right-hand sides of the above three formulae by factors of l/(s/2e) l/(V2e 2 ) and l/(2e 2 ), respectively, 
and Long's similarity variable is f/(v / 2e) times our definition of <f>. Since these factors of \/2 and e are 
completely unnecessary we eschew them in our analysis, except for the purposes of comparison with previous 
results which make use of Long's definitions. In terms of our more-orthodox definition of the small parameter 
e, the key variables in our formulation (in terms of independent variable 4>) and Long's dependent variables 
/, <?, s and independent variable y are related simply as 

/(<£) = /(y), g{4>) = -g{y), p(0 = --*(»). (8) 

e e 

where s(y) is the scaled, dimensionless static-pressure-related variable in Long's notation. 

Rewriting the three momentum equations in terms of the new variables, and also differentiating the 
z-momentum equation (|3dl) once and using the -/"-momentum equation (|3b[) to eliminate the pressure term 
from the result yields the following system: 



4> 3 p' - g 2 =0, (9a) 
4>g" + {f-l)g' + (n-l)f'g = Q, (9b) 
3 /"" + (/ - 2)^ 2 /'" + [(4n-W 3(f-l)](4,f"-f')+2[n<j>g , + (n-l)g]ct>g = 0, (9c) 

where, for the sake of subsequent analysis, we have also rewritten things entirely in terms of n, the viscous 
layer growth-rate parameter. Note that for n — 2, the /, g form of the ^-momentum equation (|9b[) can be 
integrated once: 

4>g' + (/ - 2)g = ke, (10) 

where kg is the integration constant. However, this proves not terribly useful since in practice n = 2 
appears to be far beyond the range of n in which the resulting flow solutions are physically reasonable. 
The pressure-eliminated z-momentum equation (|9cj) can be exactly integrated for two particular values of 
n, namely ^ and 1, which can be seen by recasting (|9c|) in the following form: 



it/ (»-2), T 4(n-l)(2n-l) ,, 
<^*' + - = — -/' 



I") 



where VP is defined via the rather unenlightening (and ultimately unimportant, except for this exercise 
in evaluation of integrability) relation 

* = ff" + ff>) + (13n 2 -14n + 4) + „ + _ l? + a> 

n 

The only important aspect of the definition of ^ is that it makes it easy to see that for n = i and n = 1 
the right-hand side of equation (fTTj) vanishes identically, allowing the left-hand side to be integrated. For 
n = \ (Blasius-like viscous layer growth, zero axial pressure gradient) integration of equation (|9"c)) yields 

2 /"' + (/ - 1) W" - /') + \t>9 2 = Ccb 3 - 0, (13) 

where the integration constant C can be shown to be zero by direct substitution of the definitions (J7J , 
along with the setting of n — i, (m = 0) into the primitive- variable z-momentum equation (|3d[) . 
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For n = 1 (conical self-similarity, Long's case) integration yields 



4> 2 f" + (/ - 3 W + [<Pf - 3(/ - 1)] /' + <f>g 2 = D4? = 0, (14) 

where the integration constant D can again be shown to be zero, although the analysis required to do so 
is bit less trivial than in the case n = i (cf. [Lee(1998)| .) 

The combination of equation (|9bj) together with either (|9c|) . (fl~3|) or ([T4|) constitutes a coupled system 
of two equations in two unknowns which can be solved as it stands, but to complete the comparison with 
Long's formulation we integrate (TT4"]) once more, using that g 2 — <\?p' to obtain 

M" + (/ - 1)/' + 3 P = W = 0, (15) 

where the integration constant E again must be zero due to the requirement that the nondimensional 
pressure p vanish at infinity. We have thus recovered Long's equations. 

In the general case, equations (|9al c) constitute a seventh-order system, fourth-order in /, second-order in 
g and first-order in p. For the case of the generalized Long-type solutions, i.e. on a semi-infinite domain, we 
require all velocity components to vanish at infinity - as a consequence of this, the nondimensional pressure 
must vanish at infinity as well, except in the special case n = h, where its far-field value is an arbitrary 
constant (hence can also be chosen zero with no loss of generality.) Thus, in terms of the dependent variables 
/, g and p, suitable boundary conditions are 

<fi = :/ = /' = g = 0,/" = w axls ; (f) -t oo : / ~ f as ,g ~ g as ,p -> 0, (16) 

where f as and g as refer to the asymptotic form of these variables in the far field, which must be known (at 
least to leading order) for our high-accuracy numerical solutions method to work; the far-field asymptotics 
are the subject of the next section. 

However, before getting into the details of the asymptotics and numerics, we first briefly discuss why 
Long's solutions fail the total-pressure criterion for physicality by way of comparison with the generalized 
Hall-type flows studied by Mayer and Powell, and thus further motivate an analogous generalization of the 
Long-type solutions. We first define the static and total pressure coefficients in a slightly generalized version 
of the typical fashion, which applies equally well to both finite and infinite-domain flows: 

_ P~Pe r _ PO ~ POe n „s 

p - \ P wr po \ P wi ■ [ n 

Here the subscript e denotes some chosen "outer edge" location, which can be finite but need mainly be 
sufficiently large that the vorticity of the flow solution there be negligible compared to its maximum value 
in the inner viscous core. In figure [1] typical distributions of C Pa versus or both the generalized Hall-type 
flows studied by Mayer and Powell and for Long's infinite-domain conical solutions are shown, all for flows of 
similar strength as measured by the vortex-axis value of axial velocity excess. For the generalized Hall-type 
flows all the C Po distributions are asymptotically constant in the far field and decrease monotonically as 
4> — > 0, but for Long's solutions, C Po achieves a maximum on the vortex axis, as explained previously. As 
is clear from the similarity equations, the value of the pressure (and not just its derivative) plays a role for 
nonzero axial- flow exponents m (n ^ i). For the Hall- type solutions one can set p independently of the 
velocities at the finite outer boundary <p e , but for Long- type solutions one can only specify that p vanish as 
4> — > oo and then take whatever total-pressure distributions one gets from satisfying the similarity equations 
for the given on-axis value of axial velocity; it is clear the figure (and from analogous plots of other Long- 
type flows for all admissible values of flow force, not shown here) that at least for n = 1 what one gets is a 
nonphysical C Po profile. 



4 Far-Field Asymptotics 

In order to understand the far-field behavior of the generalized Long-type flows and - as it turns out - to 
permit application of the spectral numerical scheme described in section <j5]to solve the similarity equations, 
the far-field asymptotic form of the stream function / and circulation g for general values of n is required. 
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For n = 1, following Long and including our slightly different but equivalent definition of the independent 
variable 4> one finds that 



/ ~ -Zj= + a + a^f)- 1 + (18a) 
g ~ 1 + exponentially decaying terms. (18b) 

(A detailed proof that the non-constant terms in the far-field expression for g do in fact vanish exponentially 
fast for the n = 1 case can be found in [Lee(1998)] . Accordingly, all velocity components decay as C((/) _1 ), 
and the pressure as 0((f>~ 2 ). The coefficient a\ is a continuous function of axis value of axial velocity w ax i s , 
but begins to drop with ever-increasing rapidity as this parameter becomes slightly less than zero, and it 
appears there are no solution below a limiting lower-bound values of w ax is which lies somewhere between 
-0.05 and -0.10, and coincides with a negative divergence of the ai coefficient. 

In the general case, perhaps the simplest way to find the leading-order power of <j) in the asymptotic 
series for / is by requiring the radial velocity u to decay at infinity (a property of the n = 1 solutions, which 
we also require in the general case), i.e. that 

f 

u = nf - '— -> as 4> -> oo, (19) 

and / ~ <fin is the largest power for which this is true. (For a more detailed analysis as to why this value is not 
merely an upper bound on the leading-order asymptotic-series-term exponent but in fact the unique proper 
choice, based on analysis of the inviscid terms in the scaled self-similar momentum equations, cf. |Lee(1998)] ). 
The relation between the leading-order terms in the / and ^-expansions can be found by integrating the 
inviscid form of the ^-momentum equation (|9b[) to yield g oc f 1 ~ n ~ With some further work (we 

again refer the interested reader to the thesis of Lee for the details here) , it can be shown that the proper 
generalized form of the far-field asymptotic series for / and g are ones in which each successive term is a 
power <ff~ « smaller than the preceding one: 

/ ~ fas = a_i0" + a + ai0~" + a 2 </>~™ + (20a) 

g ~ g as = k-i^- 1 ' + b^- 1 + M ( ~ + M ( ~ + (20b) 

Before developing expressions relating the undetermined coefficients in these series, recall that for n — \ and 
n = 1 the z-momentum equation (|9cl) is intcgrable. This raises the question: which form of the z-momentum 
equation should one use in the asymptotic analysis? This question turns out not have the trivial answer of 
"either" . For n — \ the above generalized far-field asymptotic series take the form 

/ ~ a _!0 2 + l + ..., g~ </, + .... (21) 

The corresponding swirl and axial velocities have the form 

« = |~1, w = 4~2a_!, (22) 

<P 9 

and thus do not decay at infinity. (For n > \ they in fact diverge.) Thus, it appears that solutions with 
vanishing far-field velocities only exist if n > i, that is, if the external axial pressure gradient is adverse 
(m < 0), and from here on we shall consider only such cases in the context of Long- type solutions which 
are potentially physically realizable, although we can find and examine solutions with nondecaying velocities 
using our hybrid numerical/asymptotic solution procedure, as well. 

For n = 1 one should use the integrated third-order form (|14[) of the equation. This is not merely for the 
convenience of using a lower-order form of the equation; as it happens, if one tries to use the fourth-order 
form of the equation when n = 1, one cannot solve for the leading-order coefficient in the /-series, since 
the relevant terms in the asymptotic hierarchy cancel identically. When one takes n — 1, substitutes the 
expansions for / and g into equations (|9b[) and (I14[) and collects terms at various orders in <p, one obtains 
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the following expression hierarchy for the coefficients, which can be solved one line at a time: 



1 - 2a 2 _ 1 = 0; 
a_i& = 0, 
2a_i6i - 6 (3 - ao) = 0, 
3fl-i&2 — 2&i(4 — a ) + aibo = 0, 
4a_i6 3 - 3& 2 (5 - a ) + 2a\b\ + a 2 b Q = 0, 



3(a - l)a_i - 2& = 0; 
&o + 2&i = 0; 

5 [ax (3 - oo) - a_ia 2 ] - 2(6 &i + & 2 ) = 

?] - 2(6 6 2 + 



(23) 



6[2a 2 (4 



ao) 



2a_ia3 



63) - b\ = 0; 



from which we see immediately that a_i = \/\J~2, ao = 1, ai is indeterminate, a 2 = 2\/2ai, 03 = ai(12 — 
ai/v2), and similarly, all higher-order a-series terms are dependent on the a priori unknown a\ term - in 
terms of our soon-to-be-described solution procedure, a\ is the eigenvalue which must be found by matching to 
the details of the inner core solution. All the b-coefficients are zero, another indication of the aforementioned 
exponentially-fast asymptotic approach of circulation to its value at infinity. 

In the more-general nonintegrable case one encounters a similar scenario, but now a_i is the undetermined 
coefficient and all higher-order coefficients are functions of it. After some rearrangement the first few /-series 
coefficients are: 



«o 



«2 



2n z {\ - n)(l - 3n) - a 2 _x(l - 2n) 2 (l - An) 
a 2 _ 1 n(l - 2n)(l - in) - 2n 4 (l - n) ' 

b n 2 [n(a a - 3) - a~i&o] 
2(l-n)(4a 2 _ 1 +n 2 ) ' 
N 2 i - N 22 



(24) 



where A^i 



30^.^3 - 2n)(l + 4n) + 6n 3 (l - n) ' ' 



a_iai(l + 2n)(l +4n) 2 - a + - 



2n J 



3a_i6 foi + (1 + n) I 3 — ao 



- 61 + oi6q(1 
n ' 



2n) 



Note that all higher-order terms in the two generalized asymptotic series can be expressed strictly in terms 
of a_i, but as this quickly leads to extremely unwieldy expressions, we have instead chosen to write each 
higher-order coefficient in terms of convenient combinations of lower-order terms in order to keep the notation 
reasonable. This kind of recursive cascade of coefficients also proves convenient from the standpoint of 
numerical implementation - the dependencies among the terms as written merely imply that they must be 
calculated in a particular order, i.e. provided with a suitable value (or initial guess) for the leading-order term 
a_i, we calculate in order ao, bo, a±, 61, and so forth when we update the coefficients after each iteration. 
The first few coefficients in the generalized far-field g-series are 



b-x 
bo 

bx 

b 2 



= 1, 



(1 - n) 



1 

2a_i 
1 

3a_i 



ao H 3 

n 



n(3 — ao)bo + 2ax(l — ri) 



axb {l - 2n) - (1 + n) ( a H 3 ) 61 + 3a 2 (l - n) 



(25) 



This difference in indeterminate coefficient seems strikingly at odds with the integrable n — 1 case, especially 
since as n varies we expect a continuum of solutions. Reassuringly, in the general case the dependence of 
a_i on the details of the viscous core becomes less and less sensitive as n — > 1 until at n = 1, numerical 
solutions of the full (nonintegrated) equations for any value of vortex strength (as measured by either flow 
force or w ax is) all yield a single numerical value for a_i = 1/v2j which is exactly the unique value found 
from doing the asymptotics using the integrated form of the z-momcntum equation for this special case. 
This behavior is illustrated in figure [21 Thus, even though the undetermined coefficients appear in different 
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Behavior of the coefficient a_i in the asymptotic series fas with respect to n for various values of 



places depending on whether n = 1 or not, the two forms agree for the unique value of n for which both 
are applicable. These formulae show that when n = 1, except for a_i not being determinable a priori, the 
coefficients given by the general-n asymptotics agree with those for the n = 1 integrated-equation results if 
one sets a_i = 1/V2. Note also that the factor (1 — n) in the numerator of bo and the denominator of a\ 
makes a\ indeterminate when n = 1. Thus the above are only strictly valid for 1. However, as we shall 
see, the above coefficients are consistent with those for the n = 1 case, in that the (in general indeterminate) 
coefficient a_i converges to a_i = 1/a/2 for all inner-core solutions (i.e. vortex strengths) as n — > 1. 

Lastly, before describing the numerical method used by us to solve the generalized equations, we note that 
there is one additional "asymptotic surprise" in store. Inspection of the expression for ao in ([24)) indicates 
that it is possible for the denominator to vanish for values of the growth-rate parameter | < n < 1 if the 
value of the leading-order coefficient a_i — (which depends continuously on the vortex strength, as measured 
by the value of u>(0)) happens to take on the value 



2n 3 (l 



(l-2n)(l-4n) 



(21) 



where the argument of the square root descends from a positive divergence as n increases from 0.5, and is real 
and positive for all n lying strictly between one-half and unity. For n = 1 it might appear at first glance that 
ai similarly blows up, but substitution of the expression for b reveals that the zero denominator is canceled 
by an equal-order vanishing numerator, leaving a finite result, albeit an indeterminate one. Note also that 
for the above values of a_i, the numerator in the expression for a takes the value — 2n (1 — n), so the zero 
denominator is not balanced by a vanishing numerator unless n = (outside the range of consideration) 
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Figure 3 - The singular behavior of ao for w aX is = 1.0, 0.5, 0.2, and —0.07. 

or n = 1, the conical Long's case. Some of this singular behavior of ao, here plotted for constant vortex 
strength and with varying n, is illustrated in figure [3j here by plotting the evolution of ao for various fixed 
values of w ax is as the vortex-growth-rate parameter n is varied. (The line segments connecting the positive 
and negative-divergence portions of each trend line are artifacts of the plotting software, but useful ones in 
terms of connecting the two distinct asymptotic branches of each curve). We were in fact first alerted to this 
possibility by way of seemingly curious behavior of the numerical iteration procedure for values of vortex 
strength which happened to correspond to values of a_i close to the above "blowup value" and which thus 
resulted in near-singular behavior in ao and the other terms depending on it. Interestingly, the occurrence of 
such "coefficients crises" appears to be closely correlated with the appearance of near-axis overshoots of the 
total pressure coefficient of the solutions, that is, with incipient nonphysical behavior. It is possible that the 
coefficient singularities define separatrices dividing physically possible from nonphysical flow regimes, with 
the "physical" portion of the solution space lying on the lower- n side of each such divergence point. 

4.1 Conservation of Axial and Angular Momentum Flux In a Constant-z Plane 

For any flow field of the type under consideration here, since there is no axial external impulse or torque 
acting on the flow, one might at first blush expect that the axial and angular momentum fluxes through any 
constant-z plane should be constant for the flow to be be physically plausible. This is indeed the case for 
Long's conical solutions, but as it turns out, does not hold for the general nonconical similarity solutions. 
However, this appears to be a case in which a superficially reasonable-seeming criterion indeed proves not so, 
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and the fact that the above properties hold for one special value of the core growth rate parameter n (which 
is also exceptional in numerous other ways, as we have detailed) turns out to be a quirk of the asymptotic 
solution properties which obtain in that case, rather than anything having a bearing on physicality of the 
solutions. Because the analysis involved here illustrates some further important aspects of the flows and the 
similarity solutions, we shall devote a page or two to it. 

To obtain the axial momentum flux through a constant- z plane, we integrate the z-momentum equa- 
tion (|3d|) with respect to the radial and polar ((f), 6) coordinates of the similarity formulation: 



[uw' + w(mw — ntfiw' )}(f)d(f) 



On using the continuity and r-momentum equations equations 
rangement, we obtain: 



— — + nepp — 2mp 



(22) 



a[b), integration by parts and some rear- 



2(m + n) / (p + w 2 )(j)d(j) = \<pw' + n<p 2 p — <pw(u — n<pw)\ 
Jo 



(23) 



Since the right-hand side of this last expression is zero at = 0, we need only consider the limit at infinity, 
for which we can make use of the asymptotic series for stream function and circulation for general n which 
we derived in the preceding section. Using that u — nf — f /4>, v — g/<f> and w = f /4>, we have: 



2 (to + n) I (p + w' 
Jo 



2(1 - In) 



-2) 



(24) 



The first thing to note about this expression is that is shows why the far- field asymptotic constant a_i 
tends to a fixed value as n —¥ 1 (and thus m — > — 1: For (to + n) —¥ 0, the entire left-hand-side of the 
above expression vanishes, hence a_i must tend to l/y/2 (the value which makes the leading-order term of 
the right-hand side vanish) , independently of whatever core- flow solution details are contained in the terms 
appearing in the left-hand-side integrand. The above integral is - up to a constant involving the common 
powers of e and z which we eliminated from both sides of the scaled z-momentum equation during the 
similarity formulation - equivalent to the axial derivative of flow force. It is clear that the leading term on 
the right-hand side of the above asymptotic expression only decays at infinity for n > 1. For n — 1 the 
power of (j) of this term is precisely zero, but the coefficient vanishes since m + n does. For | < n < 1, on 
the other hand, this term diverges unless 



0-1 



2(2n- 1) 



(25) 



The curve described by this equation in fact passes through the physically-plausible region of the solution 
space based on the C Po criterion. Thus, if we insist on constant axial momentum flux the best we can say is 
that for n < 1 there are at least some points in the physical region (based on total pressure coefficient) which 
also satisfy the requirement of constant axial momentum flux. But it is in fact not reasonable to impose 
such a requirement, and further it makes no sense that there be nonphysicality related to conservation of 
momentum, because the similarity formulation by construction conserves momentum, at least up to the order 
of the neglected terms in the small parameter. Regarding the artificiality of the constant-axial-momentum- 
flux criterion, we quote Morton(1969)| (emphasis his): 



We note that the pressure Poa — po on the axis of a rotating core embedded in an externally 
still environment at uniform pressure poo is constant only in a cylindrical field v — v(r) of 
azimuthal velocity, and that there will always be axial pressure gradients (with corresponding 
axial velocities) in rotating cores exhibiting axial development of the azimuthal velocity field 
v = v(r,z). Thus axial pressure gradients are the rule in all swirling flows suffering progressive 
lateral spread by either viscous or turbulent diffusion ... In strongly rotating cores the pressure 
thrust is of comparable importance to the momentum flux, and the axial pressure gradients 
normally formed result in a progressive transfer of force between the momentum flux and pressure 
fields ... in the absence of body forces the axial changes in flow force are due to entrainment of 
ambient axial momentum and to outer viscous stresses. 
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In the case of variation of angular momentum flux in the axial direction, we obtain an analogous result - 
the flux is constant for the special conical case due to the exponential approach of the circulation to a 
constant in the far field and corresponding vanishing of the higher-order 6-coefficients in (I20b|) , nonconstant 
in the general case, but satisfying the momentum conservation equations up to the neglected terms in the 
similarity equations. As to the physical mechanisms for the nonconstant angular momentum flux, again 
quoting Morton: 

...the axial rate of change in the flux of angular momentum (due primarily to convection, 
and in negligible proportion ... to viscous diffusion) is the result of an effective torque due to 
entrainment of azimuthal momentum at infinity and to the moment of viscous stress at infinity. 



5 Numerical Method 

We desire our numerical solutions to be of sufficient accuracy to serve as base-flow profiles for a subsequent 
numerical linear stability analysis, thus we need the basic-flow solutions to be accurate to roughly the level 
of round-off error in the latter - typically that means to the level of IEEE 64-bit floating-point arithmetic. 
(In other words, a standard kind of finite-difference approximation scheme simply will not do). To that end, 
the general similarity equations (|9b[ c) are discretized using a Chebyshev polynomial basis and a standard 
interior-plus-endpoints meshing (cf. |Boyd(1989)| ) of the Chebyshev interval x 6 [—1,1], along with a 
suitable mapping transformation to the semi-infinite "physical" (that is, <p) space: 



J 

71- 



N 



and (pj = A \ + X J , for j = 1, . . . , N, (26) 



where N is the number of spectral collocation points and A is a computational parameter (typically order of 
unity) we dub the "half-points radius" because it coincides with the physical location inside which precisely 
half the collocation points are located. The parameter A can be tuned to give a desirable level of resolution 
of the immediate near-axis region, according the observed behavior of the solutions for some intermediate 
level of numerical resolution. We use the rapidity of convergence of key numerical solution parameters to 
guide the choice of A - the speed of convergence of the asymptotic-series coefficients a±\ with increasing 
N, for example, appears to be a good guide to the choice of an optimal value of A for a given set of flow 
parameters. 

This nonlinear one-parameter mapping is designed to spread the outermost collocation points (which 
have a cosine-type clustering in computational space, i.e. are most densely clustered around the endpoints 
x = ±1), far apart as 4> — > oo, while still maintaining good spatial resolution near the rotational-symmetry 
axis of the flow, that is the origin of our radially semi-infinite coordinate system. Our preliminary numerical 
experiments showed that this stretched-grid method will only converge if the dependent variables being 
solved for decay at infinity. A little analysis shows that this is a direct consequence of the above nonlinear 
mapping between x and (/>-space being singular at x = 1, as a result of which the dependent variables 
must thus decay as <f> — > oo sufficiently rapidly to cancel out this singularity in the mapping function. The 
generalized asymptotic series for / and g thus prove not only useful but in fact necessary for our full-field 
computations. We rewrite equations (|9bl c) (substituting equation (fT5|) for ([9"c]) in the special case n = 1) in 
terms of modified functions F and G, defined as 

F{4>) := f{4>) - f as (<f>)s(<f>), G{4>) := g{4>) - g as (<f>)s(<f>), (27) 

where s(4>) is a suitable interpolation function, designed to transition smoothly from zero at 4> = to unity as 
4> — > oo, with exponential-asymptotic behavior in both limits. This gives the desired decay of the dependent 
variables at infinity, while still allowing one to use the original inner boundary conditions on / and <?; that 
is, the boundary conditions in term of the modified dependent variables are 

= : F = F' = F" - w axis = G = 0; F, G — > as — > oo. (28) 

We find the following interpolation function to work well in practice: 

. . 1 + tanh[z(x)l , Cx . . 

s((p) = , where z = —^=^= maps [—1,1] to (—00,00). (29) 

2 VI — x 2 
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Here, C is an order-unity constant which controls how quickly the function s transitions from near-zero 
values to near-unity values. We want s to approach unity exponentially fast as <f> — > oo, but we do not 
want the transition region about the inflection point to be too narrow, in order to maintain good numerical 
resolution there. After transforming the original system of equations for / and g into the modified equations 
for F and G, we set up a Newton-type iterative scheme by linearizing the variables we wish to solve for, 
i.e. the modified stream function and circulation F and G and the asymptotic coefficient a±\ (the ± sign 
on the subscript depending on whether n = 1 or not, respectively) about their values at the hypothetical 
current (fcth) iteration: 

F (k+i) =F (k) +§F (k)^ G^=G^+5G^ k \ =a%l+5a%l (30) 

where an iteration cycle begins with some initial guesses for the functions being solved for F^ and G*- * 1 
and for the numerical parameter a+|. Substitution into the modified equations and neglect of quadratic 
J-terms leads to system of linear differential equations for the functions 6F and 5G (where we drop the 
iteration superscripts to simplify the notation), with 5a±\ appearing as an eigenvalue and the fcth- iteration 
residuals (which are just the modified equations written entirely in terms of the known current-iteration 
values FW, G^ and a^) appearing as forcing terms. Next, we expand each of the two unknown functions 
in an TV-term-truncated Chebyshev expansion: 

N-l N-l 

5F{4,) = J2 CiTMti)' SG W = E d i T M<f>)), (31) 

3=0 j=0 

where Tj is the jth Chebyshev polynomial, which one can define conveniently via the trigonometric identity 
7} (cosy) = cos(jy), although in practice we use the standard 3-term recurrence 

Tj + i(x) — 2xTj(x) — Tj-i(x) for j > 1, with Tq — I and T± — x, (32) 

in order to generate the numerical values of the Tj and the needed derivatives thereof at the collocation 
points. Note that if one had reason to believe that there were a multiplicity of solutions at given flow 
parameters (with solutions as parameterized here, i.e. with the non-uniqueness resulting from use of flow 
force as the key parameter removed) it would be worth considering a more-general eigensystem approach, 
but ones such as are used to investigate multiple eigenmodes in linear stability analyses seem inapplicable 
due to the nonlinearity of the equations in question. However, as we have found no evidence of multimodality 
in our Wams-based parameterization (e.g. via any of the quite different kinds of initial guesses we tried early 
in our investigations at the same values the flow parameters w ax is and n converging to different solutions) , 
this iterative single- value technique for finding both a±i and the Chebyshev-expansion coefficients of the 
associated functions F and G proved sufficient for the task at hand, in addition to being computationally 
fast. 

Requiring the expansions to satisfy the linearized equations at each of the N — 2 interior collocation 
points, along with the five boundary conditions (|28p . also rewritten in linearized <5-form, yields a total of 
2 x (JV — 2) + 5 = 2N + 1 linear-algebraic equations for the 2N + 1 unknowns (the 27V expansion coefficients 
Cj and dj and the iterative correction to the asymptotic coefficient 6a±±), hence a well-posed linear system, 
which is solved by standard numerical techniques - we use Gaussian elimination with full (row and column) 
pivoting. We perform the computation in 128-bit emulated floating-point arithmetic with roughly 34 decimal 
digits of precision, and stop the iteration when all terms in the discrete residuals vector are less than 10 -25 
in magnitude. Beginning with N = 60, we repeat the entire iteration sequence (assuming it converges), 
each time with a larger N (increasing by 20 up to 120, then proceeding in increments of 50 from N = 150), 
using the converged numerical solution from the next-smaller N as an initial guess, until the values of the 
asymptotic coefficient a±i obtained using successive values of N agree to at least 10 significant figures. 

Convergence histories of a\ for three sets of numerical trials, all for the same flow parameters (n = 1, 
w ax is = 1) i- e - Long's vortex with a central axial velocity of unity, and a collocation mesh with half-points 
radius A = 5) are shown in Table [TJ The three trials differ only in the number of leading terms of the full 
asymptotic far-field series for / used to construct f as . The data show that in each case, the value of a± 
tends to the same constant as the numerical resolution is increased, and that retaining higher-order terms 
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-0.573422092266259 
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250 
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2.0c-10 


-0.573422092321949 


4.1c-15 


-0.573422092321946679 


7.6e-18 



Table 1 - Convergence of the asymptotic coefficient a\ with increasing spectral resolution, with respect to the 
high-resolution reference value (ai) rof = -0.57342209232194667467 . . .. 

in the far-field asymptotic expansion (which has as its price an increase of the algebraic complexity of the 
modified equations and their discretization - we refer the reader to the source code listing in the appendix of 
|Lee(1998)] , especially the segment of the listing beginning on page 142, for the gory details) generally leads 
to faster convergence, especially at the higher values of N. (Note however the slightly anomalous behavior of 
the 0(1/ (f> 3 ) approximation in the range N = 150 to 250, for which the respective relative errors are less than 
for the 0(l/4> 4 ) asymptotic series truncation, for reasons unclear.) In the table, the number of significant 
figures was computed by comparison with the reference value a\ — —0.57342209232194667467... yielded 
by a trio of ultra-high-resolution runs with TV = 400, 450 and 500 and asymptotic /-term up to 0(l/</> 4 ), 
which agreed to at least 20 significant figures. For the purpose of our general study of these solution families, 
however, 5 to 6 significant figures in the requisite asymptotic constant a±\ appears to be more than adequate, 
so for our general-n runs we have retained terms up to 0(4>~~) in f as and 0(4>^~^ 1 ) in <? as , that is, the 
leading four terms in each of the generalized far-field asymptotic series (|20b[) . 

On a typical gigahertz-class workstation with reasonably efficient emulation of 128-bit floating-point 
arithmetic a full convergence cycle (i.e. iterate at fixed iV until convergence of the numerical residuals is 
achieved, increase N and repeat until a±\ converges to at least 6 figures) for cases similar to those described 
above, with spectral resolutions ranging from N — 60 to 300, takes under a minute. The higher-resolution 
runs are where most of the expense is incurred (since the work in solving the linear system at each iteration 
scales as 0(N 3 )), but the extra work per iteration is partly offset by the fact that the higher-resolution 
runs, being bootstrapped from ever-higher-quality converged solutions for the next-smaller N, need just 
one or two iterations to converge. Low-resolution cases may take several tens of iterations, or even a few 
hundred in generalized Long's vortex cases where a good initial guess is not available or one is near the 
asymptotic-coefficients singularities mentioned above and the iterations need to be severely under-relaxed 
(e.g. by updating as F^ k+1 ^ = +ujSF < ' k \ where the relaxation factor is in the range < uj < 1) in order 
to achieve convergence. 

6 Discussion of Solutions of the Generalized Long- Type Flows 
6.1 The Solutions for n = 1 

We first examine the solutions for Long's vortex (n = 1), first for the special case of a vortex with a core 
axial velocity of unity. Unless otherwise noted, all of the plots represent solutions calculated with 200 basis 
functions for each of F and G and terms up to 0(1/ 4> 3 ) retained in the /-series, which imply solutions 
converged to well within "plotting accuracy". The modified and unmodified stream functions F and their 
first two derivatives are plotted in figure[4] As expected, F and F' are both zero at the axis, F"(0) is equal to 
the desired core axial velocity, and F and all its derivatives vanish at infinity. The same qualitative behavior 
holds if we retain more or fewer terms in f as , the only difference being that the modified functions F decays 
more or less quickly in the far field. Analogous plots of the modified and unmodified circulation appear in 
figure 

The distributions of the various velocity components for solutions with various values of core axial velocity 
ranging from —0.05 to 1 are shown in figures [6][8l by way of our boundary-condition parameterization we 
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Figure 6 - Axial velocity distributions of Long's flow (n = 1) for various values of Wans 



control the axis value of axial velocity but not whether it is a local maximum or minimum, i.e. whether the 
axial flow is jet-like or wake-like. The axial velocity (fig. [6]) is clearly jet-like for w ax i S = 1 and 0.5, begins 
to show a core axial velocity deficit for w ax is = 0.2, and for w ax is = —0.05 exhibits a region of reversed flow 
in the core. 

The azimuthal velocity (fig. [7]) shows the expected transition from a solid rotation at the center to a 
potential-vortex-like decay in the far field, with a peak value that becomes smaller as the core axial velocity 
is reduced, but which remains positive even if w a xis is less than zero. The swirl velocity distributions are 
much less sensitive to the flow parameters in the far field, in fact they all asymptote to the same curve, 
which is strictly determined by the normalization of nondimensionalized circulation in the far-field. 

Owing to the requirement of conservation of mass the radial velocity curves in figure [8] show qualitatively 
different behavior for the different types of flows: if the axial velocity is jet-like in the core, the radial velocity 
is of positive sign near the axis, becomes negative in an intermediate region and tends to zero from below 
in the far field. If the axial velocity is wake-like in the core, the radial velocity is negative near the axis, 
becomes positive in an intermediate region and tends to zero from above in the far field. 

As far as what happens when one tries to extend the parameter range further, this is quite different for 
the jet-like and wake-like solutions. For both, driving the core axis velocity further in the direction in which 
it is already tending (i.e. making the flow more jet-like or wake-like, respectively) corresponds to increasing 
the flow force. Increasing w a xis for the jet-like flows is not a problem, but making w ax is more negative in 
the wake-like case becomes extremely difficult for w ax i S < —0.05, and one must make successively smaller 
changes in this parameter value, and severely under-relax the iterations to have any chance at obtaining a 
converged solution. Figure 3 of Foster fc Duck(1982)| provides a partial explanation of why this is so. It 
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Figure 9 - The azimuthal vorticity of Long's flow (n = 1) for w aX is = 1.0. 



shows that whereas for the jet-like (Type I) solutions, the slope of the solution curve in the flow lorce-w ax i S 
plane appears to tend to a constant positive value, for the wake-like (Type II) solutions, dM / dw ax is — > oo 
as M — > oo, meaning that small changes in w ax i S produce ever-larger changes in the flow force (and in 
certain of the solution curves). The other difficulty with increasingly wake- like solutions is related to the 
large-flow-force asymptotic properties of Type II solutions described by (Foster & Smith(1989) : as M — >• oo, 
these solutions take the form of a ring-jet, with an ever-narrower annular region of positive axial velocity 
centered at finite radius (and circulation rapidly transitioning from a value which is effectively zero elsewhere 
to unity in this interior layer) . One can observe the beginnings of the development of this layer-like structure 
in figure [6l for w ax i S = —0.05. In conjunction with the flow force, the asymptotic coefficient ai also becomes 
increasingly sensitive to changes in the flow parameters as the flow becomes more wake- like. From a numerical 
standpoint (especially if one is using a Chebyshev-type spectral method, for which the set of basis functions 
tends to concentrate the spatial resolution near the endpoints of the computational domain), to achieve 
reasonable spatial resolution of the large-flow-force Type II solutions, one would need to devise a modified 
mapping of the numerical to the physical domain allowing one to concentrate some desired fraction of the 
collocation mesh points inside the developing interior layer. This is not tremendously difficult to do, but we 
do not pursue it further, since our qualitative conclusions regarding the nonphysicality of the conical Long's 
solutions (see the discussion of the total-pressure profiles below) applies whether the the solution is jet-like 
or wake-like, and whether the flow force is order unity or asymptotically large. 

In order to obtain at least a preliminary indication as to the hydrodynamic stability of the flow, we 
begin by observing that the circulation is a monotone increasing function of radius for all of the Long's 
vortex solutions, so by Rayleigh's criterion (cf. |Chandrasekhar(1969)|, §66) the flow is stable to inviscid 
axisymmetric disturbances. The azimuthal vorticity component (figure [9]) is everywhere nonnegative for the 
strictly jet-like flows, but for flows with a wake-like component (i.e. any region in which dw/dr = —ujg > 0) 
suffers a reversal of sign in the neighborhood of the axis, so the latter solutions satisfy the heuristic criterion 
of Brown fc Lopez(1990)| for susceptibility to vortex breakdown, at least of the axisymmetric variety. 



22 



Even with regard to the jet-like solutions, we know from the scalings in the similarity formulation that 
conical (m = —1) self-similarity implies a severely adverse axial pressure gradient, hence it would be very 
surprising indeed if such flows were stable (and we know from the extant stability literature that they are 
not.) More conclusively, we can quickly assure ourselves whether this is so (in the sense of sufficiency, 
but not necessity) via the stability criterion of [Leibovich fe Stewartson(1983)| based on high-azimuthal- 
wavenumber-mode asymptotics, which applies to generic slender axisymmetric vortex flow at high Reynolds 
number (the modes in question being subject to rapid viscous damping as Reynolds number decreases) and 
is defined in terms of our radial similarity coordinate and flow variables as being a sufficient condition for 
flow instability if the function 

G L s{<f>) = vQ'[CL'g'+ (u/) 2 ], where Sl = v/r. (33) 

is negative anywhere in the flow field. Plots of the function GLs{<t>) f° r the the n = 1 solutions (not shown, 
owing to the extremely different ranges of variation of this function for the different solutions) show it to 
always be negative in at least part of the solution domain, which guarantees that the flow is unstable to 
nonaxisymmetric disturbances as long as the Reynolds number is not small (the latter is required both to 
assure that the inviscid helical modes involved in the LS criterion are not too heavily damped, and for the 
assumption of slcnderness of the base flow in question to be justified). 

Thus, we expect the conically self-similar flows to be unstable, although they still appear to be reasonable 
from the standpoint of physical realizability (in an idealized disturbance- free world). That is, until we take 
a careful look at the pressure profiles. The radial distributions of static and total pressure coefficients C p 
and C p o , defined by (fTT)) , are shown in figures [TO] and QT] 

While the static pressure curves all show the expected minimum on the axis, the total pressure profiles 
(especially for the jet-like flows) are startling, in that they look like an inverted form of what one expects - 
instead of decreasing in monotone fashion from its far-field value and reaching a minimum on the axis (cor- 
responding to the location where the flow has suffered the most viscous dissipation), C Po instead reaches a 
maximum there. We have examined a variety of other well-known viscous vortex flows (viscous time-decaying 
potential, Burgers, Hall-Stewartson, Mayer-Powell) and observed that Long's is the only one of the lot which 
shows this type of behavior. Moreover, it is not just a particular subset of the flows which suffers from a non- 
physical total-pressure maximum in the viscous region; all of the other conical-flow solutions we have studied, 
of either Type I or II, have clearly nonphysical total-pressure profiles: for the Type I flows the total-pressure 
excess tends to occur near the vortex axis, whereas for the Type II flows it occurs away from the axis. For the 
large-flow-force asymptotic Type II flows, the overshoot is confined to the off-axis interior layer, but this dif- 
ferent spatial localization gives no reason to believe these flows are any more physical than the others. While 
this will not be welcome news to those who have invested much time and effort studying the structure and sta- 
bility of these flows, it seems difficult to believe that such a basic property of the solutions has gone overlooked 
for over forty years. But we have scoured the extant published literature on the subject and, amazingly, in 
not a single paper, from the original one by [Long(1961)|, to the later work of [Burggraf fc Foster(1977)| and 
the numerous studies of stability (cf. [Foster fc Duck(1982)[|Foster fe Smith(1989)[|Foster fe Jacqmin(1992)j 



Khorrami & Trivedi(1994) Ardalan, Draper fc Foster(1995)[ ), have we found a plot of total pressure or men- 



tion of the adverse-pressure-gradient nature of these flows. 

As it turns out, the behavior of C Po can be better understood in the context of the earlier related work 
of Mayer and Powell, whose flows satisfy the same general similarity laws and equations as those considered 
here (except for the n = ^ case where there is no axial pressure gradient), albeit with different boundary 
conditions. They found both favorable and adverse-pressure-gradient flows of a leading-edge-vortex-type 
character, but observed that solutions satisfying a monotonicity condition on the total pressure could only 
be found for favorable, neutral, and modestly adverse dp/dz, down to roughly m = —0.25 0. For m 
more negative than this, solutions of the equations could still be found, but all had near-axis maxima of 
C Po , the location of the maximum moving radially inward as the pressure gradient became more adverse, 
until eventually C Po was strictly monotone decreasing from the axis outward, much as is the case with the 
Long's vortex solutions described above. This raises the possibility that there are Long-type solutions of 
the generalized equations, representing vortical flows in less-severe pressure gradients (and thus having core 



2 The solutions for m = —0.4 plotted in [Mayer fc~P owcll(1992j] should not have been included, since they actually have 
a total-pressure overshoot somewhat away from the vortex axis, but the plot which would have shown this was inadvertently 
excluded from the relevant figure. 
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growth rates less than conical) which perhaps do have physically realistic total-pressure distributions. This 
indeed proves to be the case, and is the subject of the next section. 

6.2 Solutions of the General Similarity Equations 

One of the key numerical quantities associated with the generalized Long-type flows is the leading-order 
coefficient a_i of the asymptotic series for the stream function /. For a given set of flow parameters n and 
Waxis the value of this term - assuming a solution can be found - is determined iteratively via the numerical 
solution procedure described in section Sj5j and curves of a_i vs. n for various values of w ax is were previously 
shown in figure [2j The most striking feature of the plot is the simultaneous crossing point of all the curves 
at n = 1, where a_i = l/v2- We are able to obtain well-converged solutions for values of similarity growth- 
rate parameter ranging from 0.5 (even though the velocities do not decay at infinity for this value of n, the 
modified functions do and the value of <z_i is finite, so we have no trouble finding a numerical solution) to 
well over unity, that is, for axial pressure gradients ranging from neutral to even more adverse than those in 
Long's solutions. 

The thesis of Lee has a full variety of plots of the various velocity and vorticity components for the 
generalized flows with n ranging from just over 0.5 to greater than 1; for the sake of brevity we summarize 
only the key features here. The axial velocity plots (fig. IT2|) show a striking change as n is varied: for n near 
one-half, the axial velocity shows a pronounced retardation in the vortex core (which is what one expects 
physically for an adverse-pressure-gradient flow), but the retardation becomes less pronounced as n increases, 
and beyond an upper limiting value (which depends on w ax is, but commonly lies between 0.6 and 0.7), the 
axial velocity curves revert to the physical implausible jet-like form seen for the Type I Long's flows. (We 
illustrate using a slightly wake-like flow with w ax i S = —0.07; the qualitative trend is similar for w ax i s = 1 
except that the on-axis axial velocity transitions from a local minimum to an outright global maximum as 
n increases in that case). 

The radial velocity profiles are broadly similar, but note that the off-axis region where u is negative seen 
in Long's case is lacking for modestly adverse pressure gradients. Note that this does not imply a lack of 
physical fluid entrainment into the vortex core, which would be nonphysical as it would violate conservation 
of mass; due to the radial growth of the core with downstream distance (and the fact that u is velocity 
relative to the radial coordinate of a cylindrical polar coordinate system), the "entrainment velocity" is not 
u, but rather u — ncfyw, and this quantity is indeed everywhere negative except at the axis and at infinity, 
where it is zero. 

The swirl velocity profiles ffig. ll3[) are also qualitatively similar as n is varied, but show two major trends: 
the maximum value of v increases as n tends to one-half from above (i.e. as the pressure gradient weakens), 
and the rate of radial decay decreases. The first property means that swirl- velocity-based measures of vortex 
strength (which we do not control via boundary conditions, as we do the axial velocity) should take account 
of this n-dependent nature of maximum swirl velocity if one wishes to make quantitative comparisons of the 
flows at different n based on properties which depend on rotational intensity. 

The u, v and w velocities all reach their respective extrema at order-unity values of 4> before decaying 
further out. As expected from the asymptotic analysis, the O(cj)^ 1 ) decay rate of u is faster than the ~ 2 ) 
decay of v and w for n < 1 . 

As to the physicality of the generalized flows, total pressure coefficient distributions, again all for w ax i S = 
—0.07, appear in figure [Ml (The case n = 1 and w ax is = 1 was already shown in the discussion on physicality 
of Long's flows). For n near one-half, the profiles look quite reasonable: a pronounced deficit in the vortex 
core, and tending to a constant maximum in the far field. As n increases a small off-axis overshoot appears, 
and for larger n or larger values of w ax i S the overshoots become worse, eventually turning into the completely 
inverted-looking profiles already seen for Long's case. Note that owing to the different underlying swirling 
velocities (and thus vortex strengths, roughly speaking), it is not the absolute magnitude of total-pressure 
coefficient but the sign that is meaningful when comparing solutions with fixed w ax is and different n. 

From the cases we calculated, we list in Table [5] the range of n in which the total pressure coefficient 
satisfies the physicality criterion for various central axial velocities. We find that in general, higher values 
of w ax i s yield narrower ranges of physicality. Since the values of this parameter is specified as a boundary 
condition, such a trend implies that the jet velocity a physically reasonable flow can attain under a more 
adverse axial pressure gradient is smaller than for flows subject to less adverse pressure gradients, a physically 
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Figure 12 - Axial velocity distribution of generalized Long's type flows with w a xis = —0.07 for various values 
of n 
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Figure 13 - Azimuthal velocity distribution of generalized Long's type flows with w a xis = —0.07 for various 
values of n 
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Waxis 


2.00 


1.50 


1.0 


0.75 


0.50 


0.20 


-0.01 


-0.07 


n m ax of Physicality 


0.5425 


0.558 


0.58 


0.59 


0.60 


0.61 


0.60 


0.60 



Table 2 - The range of physicality for each case is 0.5 < n < n. 



intuitive result. The lower limit of the range of w ax i S , namely —0.07, is about the lowest w ax is for which we 
can hnd converged solutions using the previously described numerical scheme, at least without making any 
special allowance in the computational mesh for the increasingly ring-like nature of the flow. On the other 
hand, one can find converged solutions for much higher w ax i S than 2.0 easily. As the character of these high 
Waxis flows is qualitatively similar to w ax i S = 1.0 ones, we do not show them here. 

As far as the stability of the generalized Long-type flows is concerned, the preliminary (in the sense of 
a strictly temporal linear normal-mode analysis) results presented in [Lee(1998)] indicate a broadly similar 
spectrum of unstable bending-wave and nonaxisymmetric invisicid modes as has been well-documented for 
the conical case by [Foster fc Duck(1982)|, |Foster fc Smith(1989)[ , and [Ardalan, Draper fc Foster(1995)| . 
Unlike the conical case studies by [Khorrami fc Trivedi(1994)| 7~we do find some weakly unstable viscous 
modes for the generalized flow solutions. 



6.3 Relationship Between Long and Hall-type Solutions 

By way of conclusion, we shall explicitly demonstrate the relation between the generalized Long- type flows, 
the generalized Hall-type flows studied by Mayer & Powell, and the edge pressure of the flow considered in a 
finite domain [0, </> e ]. As mentioned in section [21 for flows with nonzero axial pressure gradient (to ^ 0) the 
edge pressure at a finite value of <p e has a deterministic effect on the flow solution. In fact, one can make 
a transition from a generalized Long type flow to a Hall type like flow by adjusting the edge pressure p e . 
To demonstrate this, we first calculate the solution for a case of generalized Long-type flow in an infinite 
domain, then pick a finite value of <p e and cut off solutions beyond it. We thus have the velocity and pressure 
fields within [0, 4> e ] of a Long-type flow, but finite-radius outer boundary conditions of the kind used to 
generate Hall-type solutions. Now by keeping v e and w e fixed at <j) e and adjusting p e , we can solve for the 
flow field within [0, <f) e ] in the same way [Mayer fc Powell(1992)| solved for the generalized Hall-type flows, 
and for the proper value of p e obtain the same solution as the latter authors. This transition behavior is 
most pronounced for the axial component of velocity and for the total pressure, which are shown in figures [TBI 
and 1 161 for a modestly-adverse-pressure-gradient scenario (i.e. one for which both a Long-type and Hall- type 
solution exists) having n = 0.6. The uppermost curve of each set shows the distribution of the flow quantity 
in question for the Long-type solution at the chosen arameter values. The successively-lower curves represent 
the effects of gradually lowering the finite-outer-boundary pressure value toward the value of p e = \{v\ +w^) 
of the Hall-type flows. As the plots show, as p e is lowered the axial velocity distribution transitions from the 
nearly pure jet-like profile of the Long-type solution to ones exhibiting the pronounced retardation at the 
axis already familiar from the generalized Hall-type solutions. The large set of curves in the figure shows 
the precise p e value associated with each curve; the smaller inset shows the same curves plotted all the way 
out to the chosen outer-boundary location <j) e =40. 

The total-pressure-coefficient plots show the transition from the Long-type ones having a pronounced 
overshoot near the axis - just how pronounced is clearer in the inset - to the monotone increasing profile 
of the Hall-type flow. As was summarized in Table [5J there are in fact (only just) some Long-type flows at 
n = 0.6 which possess monotone total-pressure profiles without requring such surgical intervention, but they 
exist at much smaller values of axial velocity than the cases shown here. 
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